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Abstract. Photoproduction of Df mesons has been measured in the ZEUS detector 
at HERA and compared with predictions of NLO pQCD calculations. The ratio of D*^ 
to D*^ cross sections has been compared to results from e + e~ experiments. Orbitally 
excited P-wave charm mesons have been observed in the Z?* ± 7r T final state. The 
fraction of D*^ 's originating from these mesons has been calculated and compared 
with that from e + e~ interactions. No evidence for radially excited mesons decaying to 
lO ' D* ± ir + Tr~ was found. The inelastic production of 3/ip mesons has been measured and 

| compared to LO and NLO pQCD predictions. 

o 

O : INTRODUCTION 

x . 

Charm photoproduction measurements have been performed at the HERA ep 
collider in the ZEUS detector from data taken during 1995-2000. Electrons or 
positrons with energy E e = 27.5 GeV collided with protons of energy E p = 820 GeV 
(1995-1997) or E p = 920 GeV (1998-2000). The ZEUS detector description can 
be found elsewhere [1]. 

The decay chain D s ± — > (j)n ± — > K + K~tt ± (38pb _1 integrated luminosity) was 
studied [2] as a continuation of a previous analysis of charm photoproduction [3]. 
The study of Df photoproduction provides another test of next-to-leading order 
(NLO) perturbative quantum chromodynamics (pQCD) calculations. 

Orbitally excited P-wave D mesons can decay to a D* by pion emission. Two 
of these states (_Di(2420) and -D|(2460)) have been found to decay into narrow 
states [4] with properties predicted by Heavy Quark Effective Theory (HQET) [5] 
and a third broad state has been seen by the CLEO collaboration [6]. A radial 
excitation of the D*^ with a mass of about 2.6 GeV decaying into D* ± 7r + 7r _ has 
been reported by DELPHI [7] but not seen by OPAL and CLEO [8,9]. 

Inelastic J ftp photoproduction proceeds via direct (resolved) processes, where the 
virtual photon (parton from the photon) interacts with a parton from the incoming 
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proton. In the dominant process, boson gluon fusion (BGF), the latter parton 
is a gluon. Photon diffraction to J ftp also contributes. The inelasticity variable, 
z = ^-p^-, can be used to distinguish these processes. Here P, p^ and q are the 
four- momenta of the incoming proton, J/ip and exchanged photon, respectively. 
From previous ZEUS data [10] the diffractive process dominates at z > 0.9, the 
direct photon process dominates at 0.4 < z < 0.9. The resolved photon contribution 
is expected to dominate at z < 0.2 [11]. 

Color singlet and color octet models have been used to calculate the above non- 
diffractive production processes in pQCD. For the former, the charm-anticharm pair 
(cc) from the hard process is identified with the physical J/ip state. In this model 
in leading order (LO) only the BGF diagram contributes to the direct channel. In 
the color octet model the cc pair from the hard process emits one or more soft 
gluons to evolve into the physical J/ip state. The free parameters of the model can 
be extracted from J/ip cross-section measurements and used in other inelastic J/ip 
production experiments. 

D% PHOTOPRODUCTION 

Df production was studied for: Q 2 < 1.0 Gel/ 2 , 130 < W 1P < 280Gel/,3 < 
p®" < 12GeV, \rj Ds \ < 1.5, where Q 2 is the photon virtuality, W 1P is the virtual 
photon proton center of mass energy, p^ a is the transverse momentum of the Df and 
r] Ds is the pseudorapidity of the Df. The effective mass of two opposite charge track 
combinations, assumed to be kaons, was calculated and plotted in Fig. la. A clear 
enhancement at the mass is seen. The effective mass of the combinations in this 
enhancement region and another track assumed to be a pion was then obtained. The 
peak in the Df mass region contained 339 ±48 Df mesons (Fig. lb), corresponding 
to a cross section of o~ ep -+D s x= 3.79 ± 0.59(stat)to'H(syst) ± 0.94(6r) nb. 

Distributions in p® s and rj Ds were compared with those for D*^ production [3] 
and with a fixed order NLO calculation [12] in which charm was produced by the 
BGF process. The signal is above the prediction (Fig. 2), particularly for 77 along 
the proton beam direction, as was the case for D*^ production. 

The ratio of the cross section for Df to D*^ production at HERA has been 
compared to that from e + e~ experiments, where the latter result is taken from a 
recent compilation [13] of fragmentation fractions to charm mesons (/(c — > D)). 
The results from the two types of interactions are: 

<J ep ^D s x/<Te P ^D*x = 0.41 ± 0.07lg;8i and f(c -> Df)/f{c -> D*+) = 0.43 ± 0.04. 
The strangeness suppression factor, 7 S , (the ratio of the probability to produce a 
strange quark to that to produce a non-strange quark), has also been compared 
to that from e + e~ experiments. From HERA the value of the above cross sec- 
tion ratio and the PYTHIA Monte Carlo was used and for e + e" the quantity 
2/(c — > Df)/[f(c — > D°) + f(c — > D + )] served as an estimator for 7 S . The values 
of 7 S for HERA and e+e", respectively, were 0.27 ± 0.04±g;g§ and 0.26 ± 0.03, im- 
plying consistency with universal charm fragmentation. 
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FIGURE 1. (a) M(K + K ) distribution for events inside the Df mass range, 
(1.94 < M(K + K~-k ± ) < 2.00 GeV). The solid curve is a fit to a Breit-Wigner convoluted with 
a Gaussian-shaped resonance and a background parameterization, a[M(K + K~) — 2mK\ b - (b) 
M(K+K-w ± ) distribution for events in the <j> mass range, (1.0115 < M(K + K~) < 1.0275 GeV). 
The solid curve is a fit to a Gaussian plus an exponential background. 
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FIGURE 2. Differential cross sections for the photoproduction reaction ep — > DX: (a) da/dp® 
and (b) da jdr\ D , where D stands for D* or D s . Inner (outer) error bars show statistical (statistical 
and systematic added in quadrature) errors. The D s (dots) and D* (triangles) data are compared 
with NLO predictions for D s (full curves) and D* (dashed curves) with two parameter settings: 
m c = 1.5 GeV, [1r = (thick curves) and m c — 1.2 GeV, [Ib. = 0.5m^ (thin curves). 



EXCITED CHARM MESONS 



As a basis for the study of higher excitatations of charm mesons, an enlarged 
sample of data (integrated luminosity of llOpfcr 1 ), containing a clean signal of 
D*^ mesons from both photoproduction and deep inelastic scattering, was used 
[14]. Events in the mass range 1.83 < M(Ktt) < 1.90 GeV, 0.144 < M(Knn s ) - 
M(Ktt) < 0.147 GeV were chosen (n s is the low momentum pion in the D* decay). 



The background (estimated from events in which the K and 7r in the D° mass range 
have the same charge) has been subtracted, yielding 27286 ± 232 D* ± . 

For orbital excitations an extra track, 7r 4 , was added to the D*^ candidate and 
the effective mass combination, M(Kir7r s 7r A ) - M(K7nr s ) +M(D*)(2.01 GeV), was 
evaluated. An enhancement in the mass distribution with total charge zero is seen 
in Fig. 3a. This spectrum was fitted to D® and Breit-Wigner shapes with 
masses and widths fixed [4], and convoluted with a Gaussian function with a width 
as in the Monte Carlo simulation. The background was described by x a e~^ x+/yx , 
where x = M(i^7T7r s 7r 4 ) — M(Kinr s ) — m n , with a, f3 and 7 constant. Helicity 
angle distributions for D® and Df proportional to 1 + 3 cos 2 9 and 1 — cos 2 9, 
respectively, were folded in for the fit. Here, 9 is the angle between 714 and n s in 
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FIGURE 3. Distribution of M(KTnr s n 4 ) - M(Knn s ) + M(D*)(2M0GeV) for D* and tt 4 with 
opposite charges, (a)-(b) Fit to two Breit-Wigner shapes convoluted with a Gaussian (full curve). 
The dashed histogram is for D* and 7r 4 with the same charge, (c) Fit included extra Gaussian 
with free mass and width. The dotted curves are the fitted combinatorial background. 

the D* ± rest frame. A closer look, (Fig. 3b), indicates an excess of events near 
2.4 GeV. An extra Gaussian was included to better fit the data. The fit (Fig. 
3c) yielded 526 ± 65, 203 ± 60, and 211 ± 49 entries for the number of D?, Df 
and the additional Gaussian combinations. The mass of the extra Gaussian was 



2398.1 ± 2.1±o|MeV with a width consistent with the detector resolution. 

The ratios of the rates of the D* ±r K T decay channel of D® and Df to D*^ 1 are 
3.40 ± 0.42±g;^% and 1.37 =b 0.40±g;§|%, respectively. Extrapolating to the full 
kinematic region and using [13] along with known branching ratios [4] as well as 
isospin conservation, the result is given in table 1. This result is consistent with 
those from e + e~ experiments. 

TABLE 1. Comparison of and Df production rates. For ZEUS 
errors are statistical, systematic and extrapolation errors. For CLEO, 
OPAL and ALEPH the statistical and systematic errors have been added 
in quadrature. The DELPHI results are without systematic errors. 



Experiment 


/(c-£>?) [%} 


f(c - Df) [%} 


ZEUS Prelim. 


1.46 ±0.181^ ±0.06 


2.00 ±0.581q Is ±0.41 


CLEO [15] 


1.8 ±0.3 


1.9 ±0.3 


OPAL [16] 


2.1 ±0.8 


5.2 ±2.6 


ALEPH Prelim. [17] 


1.6 ±0.5 


4.7 ± 1.0 


DELPHI Prelim. [18] 


1.9 ±0.4 


4.7 ± 1.3 



In order to search for a radially excited D* meson, D*\ the combination 
M (i^7T7r s 7r47r 5 ) — M(Kim s ) + M(D*), where 7r 4 and n 5 are oppositely charged pions 
with p± > 0.125GeV, was fitted along with a background distribution of x a e~^ x , 
where x = M(K 7171^^11^ — M(Kixix s ) — 2m n . No peak was seen in the expected 
mass range, (2.59 — 2.67 GeV) (Fig. 4). An upper limit was obtained by fitting 
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FIGURE 4. Distribution of M(K 7r7r s 7r 4 7r 5 ) - M(Ktttt s ) + M(D*) for the D** candidates with 
the D* window within the rectangle. Inset: Dashed histogram is the Monte Carlo signal normal- 
ized to the upper limit and added to the fit interpolation (dotted curve) in the D* window. 



the background outside this range, interpolating within the range and subtracting 
this from the data in the mass range. The 95% confidence level upper limit for 
the D* ± ii + 7i^ decay relative to D*^ was found to be 2.3%. Extrapolating to the 
full kinematic range and using [13], f(c — > D* + ) • B Df > + ^ D „ +7T+7T - < 0.7% at 95% 
confidence level was obtained. The equivalent OPAL limit is 1.2% [8]. 



INELASTIC J/ip PRODUCTION 

The decay channel of J/if} for 0.4 < z < 0.9, 50 < W 1P < 180 GeV and 

Q 2 < lGeV 2 , using the 1996-1997 data, has been studied [19]. Distributions in 
z, rapidity and transverse momentum squared of the J /if} and a comparison with 
theoretical expectations are shown in Fig. 5. The z dependence of the data is not 
described in magnitude by the LO color singlet and octet model with octet matrix 
elements calculated from the CDF data [20,23]. On the other hand, the NLO color 
singlet model [21] roughly fits the spectrum for p±_^ > 1 GeV. For p±^ > 2 GeV 
the LO model with octet matrix elements from CLEO data agrees with the data 
for high z only [22,24]. Currently there is no calculation in NLO for the rapidity 
distribution. The NLO calculation agrees with the p\^ data. 
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FIGURE 5. (a) z distribution for various p±^ cuts: no cut (squares), p±^ > lGeV (circles) 
and px,/, > 2 GeV (triangles). Inner (outer) error bars are statistical (quadratic sum of statistical 
and systematic) errors. Lower pair of solid curves are a prediction of color singlet and octet model 
[20] for p±,p > 1 GeV. Separation of the curves indicates the uncertainty in the color octet matrix 
elements. Upper pair of solid curves includes a scale factor of ~ 3. Dotted curve is the color 
singlet NLO prediction for the direct photon process and p±^ > lGeV [21]. Dashed curve is 
prediction of the color singlet and octet models for p±^, > 2 GeV [22] . (b) Rapidity distribution 
for p±^ > lGeV. Solid curves as in (a). Dotted curve is the LO contribution of the direct 
photon color singlet component, (c) p\^ distribution: Solid curve is the prediction of the NLO 
calculation [21]. 
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